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ABSTRACT 

A VLF relative navigation system makes use of the fact that, at 
any given point on the earth, phase delay of a received VLF signal is 
highly stable and predictable. As the receiver is physically moved, 
phase delay changes linearly with distance from the transmitting sta- 
tion, so that by keeping track of the phase delay of the received 
signal from several VLF stations, one may keep an accurate plot of 
geographical position. 

Two problems experienced in measuring phase delay are the diurnal 
shift and long path interference. An investigation is made into a 
method of predicting the diurnal phase shift and the resultant phase 
due to simultaneous reception of long and short path signals. Also 
investigated is a receiving antenna having a cardioid shaped radia- 
tion pattern which could provide discrimination against long path 


signals. 
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1. Introduction 

Certain characteristics of the Very Low Frequency portion of the 
electromagnetic spectrum make feasible many important applications, some 
relatively new and others not so new. One striking application is the 
—— of radio waves to submerged submarines which have not even 
an antenna above the surface. Another is the transmission of messages 
to mobile units operating halfway around the world. 

The application which inspired this paper, however, is that of a 
VLF relative navigation system. Such a system, using VLF transmitting 
stations already broadcasting for other purposes, has been tested and 
found feasible for both ships and aircraft operating at long distances 
from the transmitting stations. The system makes use of the fact that, 
at any given point on the earth, phase delay of a received VLF signal is 
highly stable and predictable. As the receiver is physically moved, 
phase delay changes linearly with distance from the transmitting sta- 
tion, so that by keeping track of the phase delay of the received signal 
from several VLF stations, one may keep an accurate plot of geographical 
location.. 

This paper attempts to resolve one particular problem associated 
with the relative navigation system application. This problem is the 
result of interference noted from the reception of a VLF signal which 
has travelled along a great circle path more than halfway around the 
circumference of the earth, i.e., interference of the so-called long 
path signal with the short path signal. Section 2 is a study of as- 
pects of propagation peculiar to the VLF spectrum, while Section 3 
discusses in detail the VLF relative navigation system. Sections 4 


and 5 provide information on VLF receivers and frequency standards 


applicable to the navigation system. 

Discrimination against the long path Signal is desirable to mini- 
mize navigational errors; this may be accomplished easily aboard ship 
by the use of an antenna system having a cardioid-type radiation pattern. 
VLF antennas are discussed in Section 6, with emphasis on the use of a 
loop in conjunction with a whip antenna to provide the desired cardioid 
pattern. 

Sections 7 and 8 provide a detailed discussion of the prediction 
of attenuation and phase delay of VLF waves. These predictions are 
necessary if one is to calculate in advance the effect of the long path 
signal on the observed resultant phase delay at any given geographical 
location. Section 9 provides an example of just such a calculation. 
Observed data on long path signals could not be obtained locally, since 
there were no VLF Stations at a sufficient distance from Monterey, Cal- 
ifornia, to provide a measurable long path signal. It was thought 
that station GBR, transmitting on 16.0 Khz from Rugby, England, would 
provide this data; however, this station was off the air during the 


period available for experimental investigation (July-November 1966). 
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2. Characteristics of Very Low Frequency Radio Propagation 

The characteristics of VLF propagation which make these frequencies 
useful for long distance communication and highly accurate (to within a 
few parts in 10115 frequency comparison are a relatively low path atten- 
uation which is stable with time and a phase shift which is very constant, 
and predictable, over long periods of time. VLF also overcomes (to a 
large extent) the two major disadvantages of HF transmissions, unrelia- 
bility during ionospheric disturbances and fairly rapid and deep fading 
[iie 

Ihe major disadvantages of VIF transmissions are the very narrow 
(20-150 Hz) bandwidths available and the large transmitting installations 
that are necessary. The narrow bandwidths available are due to the small 
spectrum (10-30 kHz) and, mainly, to the high Q of VLF antenna systems. 
Since these antennas are very short with respect to a wavelength their 
input impedance consists of a large reactance and a very small radiation 
resistance. To obtain maximum efficiency, therefore, it is mandatory 
that all losses in the antenna, loading coils, and ground be kept to a 
minimum. The resulting low resistance coupled with high reactance gives 
these systems an inherent high Q and consequent narrow bandwidth. 

The high powered transmitters are required to obtain large areas 
of coverage and because of the high level of atmospheric noise present 
at VLF. This atmospheric noise is due mainly to lightning discharges 
which are centered in frequency at about 11 kHz. A typical discharge 
wlll radiate approximately 10? joules and have a 3db bandwidth of about 
12 kHz [6]. 

Theory of VLF Propagation 


Ihe two theories that have been predominantly used to explain the 
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characteristics of VLF propagation are the ray theory and the waveguide 
mode theory. It has been stated [28] that the mode and ray theories are 
related and, if all modes and all rays were taken into account with ap- 
propriate corrections for losses, both methods would yield the same re- 
sults. However, at long distances the mode theory is simpler to use 
because only one mode need be considered but many individual rays must 
be taken into consideration. At short distances many modes are present 
making the mode theory calculations complicated, but only a few rays need 
be considered since the higher order rays are incident on the ionosphere 
at much smaller angles where the reflection coefficient is small. In 
this paper the only distances which will be considered are those for 
which the mode theory calculations may be assumed valid with only one 
mode present (i.e., greater than 4000 km from the transmitter). 

From observed measurements of the electric field within the earth- 
lonosphere waveguide it has been found that the model shown in Figure 
2-1 can be used to approximate the actual physical case. The mode of 
principal interest at long distances is the TM) > shown in the left side 
of Figure 2-1, the TM, and higher order modes being small enough in 
amplitude beyond 4 Mm to be neglected. From Watt and Croghan [18], Ёог 
distances greater than 1 Mm the vertical electric field can be represen- 


ted as the sum of a very few waveguide modes: 


n= l to 3 
E, = 2. : Hoe d» 1 Mm 


Here E, x Ls the electric field strength for each particular mode (in 
2 


db relative to 1 v/m) 


E n [db, v/m] = 10 Log P. + K, - 10 Log f 


E 
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where 
Р. = power radiated into the half space above the earth (watts). 
K, = 104.3 db for modes 1,2, --- a constant which relates field 


strengths at the surface to power radiated. 


f = frequency (Hz) 
h; г> effective ionospheric height at the transmitter (m). 
3 
hi E effective ionospheric height at the receiver (m). 
? 


mode excitation modification factor, defined as the ratio 


Ant 


of power launched into the concentric spherical shell guide relative to 


that for a flat guide with perfectly conducting earth. 


a earth's radius (6.4 Mm) 


d 


path length (m) 
(a sin d/a) accounts for energy spreading in the spherical 
guide structure. 
EST - mode field modification factor for a curved earth at the 
ionospheric height at the receiver. 
20 Log Gh,n = height gain factor (db) relative to the field at 


the surface. 


L 


d.n = 1055 due to discontinuities along the path (db). 
2 


effective attenuation rate in db/Mm (different for each 
mode). 
These factors are discussed at some length by Watt and Croghan and will 


not be commented on further in this paper, with the exception of the 
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attenuation factor. 


For paths with uniform boundary surfaces, the attenuation factor is 


given by 
< = MCP dK; + > ç 
where 
M(@) = earth's magnetic field factor 


= attenuation with magnetic field 
attenuation without magnetic field 


g = direction of propagation with respect to the earth's 
magnetic field 


x; = attenuation constant of the ionosphere 


I 


Xg = attenuation constant of the ground. 

Ki and Xç are related to the reflection coefficients of the earth and 
the lonosphere which are in turn related to the conductivities of these 
two media. 

The attenuation curves shown in Figures 2-2 and 2-3 are derived 
for M(@) = 1, i.e. independent of the earth's magnetic field. Here 
£ is the conductivity gradient of the ionosphere. From these curves the 
difference in attenuation of the first two modes is apparent. At 15 kHz 
and for a 4 Mm path the smallest difference in signal strength between 
the two modes is about 13 db, justifying the earlier assumption that 
only the first mode need be considered at these distances. Figures 
2-4, 2-5, and 2-6 are shown for comparison and while there are some dif- 
ferences in the attenuation curves for n = 1, the n = 2 curves show the 
same larger attenuation rates (with respect to the n = 1 curves). Fig- 
ure 2-7, while not a set of calculated curves, shows very vividly the 


effect of the earth's magnetic field on the attenuation factor. 
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The phase velocity of a certain mode for a fixed set of conditions 
can be determined from a solution of the mode equation. This is done 
by Wait [11] and Wait & Spies [23] resulting in the curves shown in 
Figures 2-8 through 2-11. The effect on phase velocity of several fac- 
tors can be seen here. Assuming a constant frequency, it can be seen 
that the height of the ionosphere has a lrge effect on the phase velocity, 
lowering it with increasing height. This is the cause of the familiar 
diurnal shift of phase delay of a received VLF signal. At night the D 
layer of the ionosphere (the lowest layer) dissipates and the effective 
height of the ionosphere is increased by about 20 km. The resultant de- 
crease in phase velocity causes an increase in the phase delay of the 
received signal as compared with a local standard. After sunrise the 
sun's rays ionize the lower portions of the ionosphere, the D layer re- 
turns, and the process repeats itself. 

Comparing Figures 2-8 and 2-9 will show that an imperfectly re- 
flecting ionosphere (e.g., an exponential variation of ion density) will 
increase the phase velocity slightly. The effect of finite ground con- 
ductivity is seen in Figure 2-10, indicating that the phase velocity 
will be slightly lower over land (as compared to sea). From Figure 2-11 
it appears that the effect of the earth's magnetic field on phase vel- 
ocity is small in comparison to the other factors. (The effect is much 
more pronounced for the n = 2 mode; however, that mode is being consid- 
ered negligible in amplitude here). Wait and Spies have many more curves 
of phase velocity, attenuation factors, excitation factors, etc., which 
would be of interest to anyone planning a more detailed study of VLF 


propagation. 
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Experimentally Obtained Data 

Figures 2-12 through 2-14 along with the experimental curves in 
Figures 2-2, 2-3 and 2-7, are examples of the vast amount of experi- 
mental data taken on the attenuation rates and phase velocities of VLF 
radio waves. The agreement of this data with calculated curves in cer- 
tain cases (e.g., Figures 2-13 and 2-14) is remarkable and provides a 
basis for using the calculated curves to predict phase velocities and 
attenuation rates. 

Figure 2-12 shows the many factors that must be taken into account 
in predicting attenuation rates: time of day, direction of propagation, 
frequency, and ground conductivity. In the 10-30 kHz frequency range, 
the range of attenuation rates is about 4 db/Mm (minimum to maximum) 
indicating that for long distances a considerable difference in signal 
Strength could be expected depending on the path traveled. 

In this paper the experimental data will be used wherever possible 
and the calculated curves will be used as a check and to fill in the 
gaps not covered by experimental evidence (e.g., attenuation of the 
n = 2 mode). 

Summary 

The one characteristic of VLF radio transmissions of interest in 
this paper is the good phase stability of the signals over long periods 
of time and the predictability of this phase. Because of this phase 
stability frequency comparisons over long distances are possible with an 
accuracy on the order of one part in 1011. тһе amplitude and phase of 
these signals can be predicted fairly accurately at distances greater 


than one megameter from the transmitter by using the mode theory of 


propagation. The actual solution of the mode theory equations is a 
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rather complex problem; however, the solutions of Wait and others have 


been published and show a very good agreement with experimental evidence 


in some cases. 
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3. A Proposed Very Low Frequency Relative Navigation System 

The proposed very low frequency (VLF) relative navigation system 
has been outlined by Stanbrough and Keilly [16], and presented in de- 
tail by Lake [25]. The system has many advantages, providing (1) rela- 
tive simplicity; (2) a long range, accurate system with worldwide 
coverage; (3) the capability of using VLF phase-stabilized stations 
already on the air, with no requirement for additional stations; and 
(4) reliable 24-hour coverage. The only required equipment for ship- 
board use of this system are one or more phase tracking receivers and, 
preferably, a directional antenna system such as a loop used in con- 
junction with a sense antenna (whip). The latter is discussed elsewhere 
in this paper. One recorder channel to provide a continuous record of 
phase for each tracking receiver used would also be helpful. 

The extreme phase stability of VLF radio waves propagated over great 
distances, plus or minus one microsecond under normal conditions, is 
the foundation of this navigation system. An oscillator at the distant 
VLF transmitter provides the station with phase-stabilized emissions, 
while a stable oscillator aboard ship at the same frequency provides 
the local reference. This local oscillator is used in conjunction with 
a phase tracking receiver to provide an indication of relative phase of 
the VLF wave at the receiver. A high-quality quartz crystal oscillator 
would be suitable for use aboard ship and a typical unit in use at the 
U.S. Naval Postgraduate School provides a frequency drift (due to crystal 


aging) of about 2.5 parts in 101! 


per day. This drift would be deter- 
mined for the particular oscillator in use aboard ship and its frequency 


would then be corrected on a regular basis, perhaps daily, to keep it 


within some specified frequency offset. At the time of taking a fix, 
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the phase variation due to the frequency offset of the oscillator would 
then be known and a compensation could be made. 

Assuming a nominal value of 3 x 10° m/sec for the phase velocity 
of a VLF wave, the number of microseconds required for the wave to 
travel one nautical mile is found to be 6.18. In view of this, it is 
clear that for each increase in phase delay of 6.18 usec, the ship has 
moved one nautical mile farther from the transmitter. 

The VLF relative navigation system makes use of the foregoing; by 
keeping track of the phase of two or more stations as the ship moves 
away from an initial known point, an indication of the ship's geograph- 
ical position is available at all times. By programming a computer to 
make the calculations and necessary corrections, one could even obtain 
a continuous track of the ship's movement. It is, however, necessary 
to use VLF stations far enough away that mode interference is not ex- 
perienced, and only the first waveguide mode is involved (greater than 
4000 km). 

One problem complicating the navigation system is that of the 
diurnal phase variation experienced with VLF waves. Fortunately, how- 
ever, this variation is predictable and the shipboard navigator can be 
provided with taBles to correct for the — À Figure 3-1 illus- 
trates the diurnal variation, and is an actual plot of the phase of the 
13.6 kHz signal from the Omega Navigation System transmitter located at 
Haiku, Hawaii, 9 to a rubidium vapor frequency standard at the 
U.S. Naval Postgraduate School. 

A first are ime SR to a solution of the problem of diurnal 
phase change may be made by assuming a trapezoidal pattern, or trapezium, 


as shown in Figure 3-2. A comparison with Figure 3-1 shows the similar- 
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Figure 3-2. Trapezium Approximation of Diurnal Shift. 

ity of the two patterns. Time intervals AB and EF have the same values 
of relative phase and represent the period of time during which the 
entire propagation path is in sunlight. Assuming the transmitter loca- 
tion west of the receiver, point B is the time of local (receiver) sunset, 
point C indicates sunset at the transmitter, point D local sunrise, and 
point E sunrise at the transmitter. It has been found that correction 
of sunrise and sunset times to those viewed from the effective height of 
the lower boundary of the D layer of the ionosphere (about 70 km during 
the day and 90 km at night) provides better correlation with the observed 
discontinuities in the phase plot. 

The diurnal shift (At in Figure 3-2) varies with distance from the 
transmitter, since it is caused by the fact that the phase velocity of 
a VLF wave has one value during the day and a lower value at night. As 
discussed elsewhere in this paper, these phase velocities may be pre- 
dicted for given paths, allowing the calculation of the diurnal shift 
at any given geographical position. 


Other, more precise, models of the phase plot may be determined 
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(as shown later in this paper). Whichever model is used, the procedure 
to be used in navigation is to predict the deviation of phase from the 
reference (all-daylight path) value, for the particular time of the 
desired fix, and subtract from indicated phase. As mentioned earlier, 
these deviations could be prepared and published in advance. 

An electronic navigation system such as that being considered is 
known as a circular grid system, as opposed to a hyperbolic system such 
as LORAN, which makes use of phase difference between a master and a 
slave station. In the VLF relative system, navigational charts may be 
overprinted with circles representing lines of constant time delay drawn 
about each VLF station as a center. These would be constructed using 
the 6.18 usec/mile discussed earlier. (It should be noted that use of 
these charts introduces a small error since the phase velocity of a VLF 
wave is not exactly equal to the velocity of light; hence, actual phase 
delay differs slightly from 6.18 usec/mile.) 

In order to use such an overprinted chart, one would first need to 
plot the ship's position from other information, which would then be 
used as a reference for future fixes obtained by the VLF relative sys- 
tem. The phase indicated by the tracking receivers for each VLF station 
being tracked would also be noted at the time the reference fix was ob- 
tained. At the time the VLF fix is desired, the phase readings on each 
VLF station would again be noted. For those stations whose propagation 
paths are all in daylight, the phase difference could then be transferred 
directly to the chart. This would be in the form of a time difference 
for each station which would indicate the distance travelled away from 
or toward the station. A short segment of a circle for each station 


would be interpolated on the chart between those segments already over- 
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printed. The intersection of these circular segments would then be the 
ship's position. 

For those stations whose propagation paths are partially or totally 
in darkness, the correction for diurnal shift corresponding to the ship's 
approximate position and time of day would have to be applied. There 
would be a correction for the time of the reference position and one for 
the time of the new position. The corrections would be subtracted in 
both cases, since phase delay is greater at night. After the diurnal 
correction has been applied, the phase information is used as described 
in the last paragraph. 

It is important to note that, since stations being used for navi- 
gation will be at least 4000 km distant, the circles of constant phase 
will appear as straight lines in small areas. In fact, if overprinted 
charts are not available, straight-line approximations of the circular 
segments may be used if the true bearings are known from the ship to 
the transmitting stations. These could be obtained from a previously- 
prepared table. The straight-line approximations will be adequate if 
the ship's position is known within about 50 miles of its actual posi- 
tion [25]. 

A very complete example of a navigational fix obtained by this 
method, and using the trapezoidal approximation of diurnal phase shift, 
is given by Lake [25]. 

There are certain problems associated with the VLF relative navi- 
gation method, some of which are (1) cycle-slipping, (2) interference 
from the long path (around the world) signal, (3) inexactness of the 
trapezoidal model of diurnal phase change, (4) transmitter outages, 


(5) frequency offset of the transmitter master oscillator, and (6) 
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sudden ionospheric disturbances (SID's). 

Cycle-slipping is a phenomenon in which the tracking receiver indi- 
cates and tracks one or more full cycles ahead or behind the true rela- 
tive phase of the VLF wave. This often occurs near sunrise or sunset 
when relative phase is changing rapidly. Frequently the diurnal phase 
shift will be more than one cycle, so that cycle-slipping is especially 
likely to secu Cycle-slipping may also be caused by interference from 
waveguide modes higher than the first order, and by interference from the 
long path signal [27]. 

The long path signal may be thought of as a phasor, having some 
magnitude and phase relative to the short path signal. Together, the 
two signals form a resultant which constitutes the composite signal re- 


ceived by an omnidirectional antenna. Figure 3-3 shows the phasor addi- 


tion of the two signals, with Ey - |E] /09 representing the short path 
signal (taken with phase - 0? for convenience), EP = |Е /ө° the long 


path signal, and Ep = Egl Z° the resultant. The notation E is used to 


indicate a phasor quantity. Since the long path signal may travel 
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Figure 3-3.  Phasor Diagram 
hundreds or thousands of miles f&rther than the short path signal, the 
angle 6 may take on any value between 0° and 360°. It is clear, then, 
that for any value of EE | which is an appreciable fraction of [Ej , the 


resultant will vary considerably from the short path signal. Since such 
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a condition will occur under favorable circumstances, the need for a di- 
rectional antenna which will discriminate against the long path signal 
is evident. The long path signal is particularly likely to be appreciable 
near the antipode of the transmitter, as well as at other locations where 
the short path signal for any reason undergoes high attenuation concur- 
rent with low attenuation of the long path signal. 

The trapezoidal model of diurnal phase change, discussed earlier in 
this section, is a good approximation in many cases; however, there are 
also frequent instances in which it introduces excessive error. A more 
refined model of diurnal phase change, useful in all cases, is presented 
in a later section. 

Transmitter outages may occur on a scheduled basis for maintenance, 
in addition to those unscheduled outages due to equipment failure which 
invariably occur. During the interval of a transmitter outage, the 
indicated phase of a tracking receiver remains constant. Since the 
transmitter master oscillator normally continues running during the 
outage, the receiver will resume tracking on the correct cycle when VLF 
transmissions resume, if the receiver has remained geographically sta- 
tionary. If, however, the ship has moved a distance on the order of 
some multiple of a wavelength, the receiver will track the incorrect 
cycle (the same one as if it had remained stationary). The problem of 
determining the correct cycle is known as lane identification, and is 
minimized in other systems by terine lower frequency (longer wave- 
length) signals. Additional information, such as the ship's dead 
reckoning track or positional information from other means, is necessary 
to identify the correct lane, or cycle. 


Frequency offset of the transmitter master oscillator of a VLF 
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Station is occasionally made, for example to bring the transmitted fre- 
quency into coincidence with the United States Frequency Standard. These 
frequency corrections would occur infrequently and advance notice thereof 
could be given to users of a VLF relative navigation system. With such 
advance notice, the navigator would have no difficulty applying a cor- 
rection for the offset. 

Sudden ionospheric disturbances are caused by three phenomena: solar 
flares, magnetic storms, and high altitude nuclear bursts. While these 
disturbances occur infrequently, they cause phase disturbances on VLF 
circuits over wide areas of the earth, causing errors in VLF navigation 
in the case of severe disturbances. Chilton, Steele, and Crombie [19] 
have pointed out that solar flares cause a depression in the effective 
height of the ionosphere, which in turn causes increased phase velocity, 
resulting in a lower value of relative phase at the receiving point. A 
typical flare has an onset time of one to 30 minutes and requires from 
one to two or three hours to return to normal. A typical phase change 
caused by the flare might be 10 microseconds, but varies with path 


length, frequency, and amount of ionosphere height change. 
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4. VLF Receivers. 

The receiver used to obtain experimental data on phase and amplitude 
of various VLF stations for this paper was a Tracor Model 599G VLF phase 
tracking receiver. A simplified block diagram showing the essential ele- 
ments of such a tracking receiver is seen in Figure 4-l. The output of 
the block marked VLF Receiver is at an intermediate frequency of 1 kHz 
with a certain phase D p? and is one input of the phase comparator. 

The other input of the phase comparator, a 1 kHz signal of phase 4, 

is derived from the external frequency standard. The output of the phase 
comparator provides the electronic phase shifter with an error signal 
A9 , which tells the phase shifter how much to shift the phase of the 
external standard to maintain phase tracking. The output of the phase 
shifter then goes to the frequency synthesizer which provides the ref- 
erence signal ( 4. The condition for phase tracking is Ag = Aur F - 

9, = 0. 

The phase tracking servo and long integration time of the receiver 
provide the capability of tracking with an input signal-to-noise ratio 
of -50 db. The rated sensitivity is 0.01 microvolts of signal at the 
antenna terminals to provide phase tracking. Using a loop antenna made 
by Tracor especially for this receiver (Model 599-600), this corresponds 
to a signal strength of 0.3 microvolts/meter. With a 150 sec servo time 
constant (front panel switch provides time constants of 5, 15, 50, and 
150 sec) the equivalent noise bandwidth is 0.002 Hz. A front-panel 
Switch on the receiver allows selection of either a broadband RF filter 
or any one of five narrowband ones. ‘The narrowband filters provide an 
RF bandwidth of 500 Hz, while the IF bandwidth is 50 Hz. The receiver 


is tunable from 8.0 to 31.9 kHz in 0.1 kHz steps. 
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The receiver has available several outputs that proved useful in 
experimental work for this paper: (a) 100 psec and 10 psec full scale 
recorder outputs, (b) a front panel bi-directional digital counter for 
cumulative phase change, (c) a 100 kHz coherent output, and (d) an AGC- 
derived relative carrier level recorder output. The 100 psec and 10 
psec full scale recorder outputs use the receiver's internal phase com- 
parator to give a dc voltage level proportional to the phase difference 
between the external frequency standard and the received VLF signal. 
(Note: the phase comparator of the receiver does not actually operate 
at 1 kHz as indicated in the discussion of Figure 4-1. The receiver 
contains two linear phase comparators; one operates at 100 kHz to pro- 
vide the 10 psec full scale recorder output, while the other operates at 
10 kHz for the 100 psec full scale recorder output.) These recorder 
outputs were convenient and proved quite satisfactory for continuous 
monitoring of the phase of any phase-stable VLF station. Moseley Model 
5100B strip chart recorders were used for monitoring these and the 
other receiver outputs, as desired. 

The electronic phase shifter of the receiver actually operates at 
10 MHz. One cycle of phase shift then corresponds to 0.1 msec drift 
which feeds the front panel bi-directional phase counter. This digital 
counter provides a cumulative indication, so that by taking readings at 
two different times, one may determine by subtraction the exact phase 
drift that has occurred during the elapsed time. 

The 100 kHz coherent output of the receiver provides a 100 kHz 
signal from the frequency synthesizer section shown in Figure 4-1. This 
signal thus has a phase corresponding to am of the received VLF signal. 


This output proved useful as an input to an external phase comparator, 
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allowing a check on the phase of any local frequency standard (the other 
input to the comparator) with respect to the received VLF signal. Both 
the rubidium vapor frequency standard and the quartz crystal oscillator 
in use at the Naval Postgraduate School have 100 kHz outputs suitable 
for this (see Frequency Standards section). 

The relative carrier level output of the tracking receiver provides 
a de level proportional to the AGC voltage developed in the receiver and 
is thus a relative indication of received signal strength of the VLF 
Station. The AGC has a range of approximately 40 db, and peaks at about 
40 db above the minimum signal necessary for tracking. 

One problem occasionally encountered in experimental work was the 
very long integration time of the receiver. For example, a cardioid 
antenna system must be adjusted for proper phase and amplitude of the 
whip antenna when changing to a new frequency. To get a good indication 
of phase and amplitude, it was usually necessary to wait five to 10 min- 
utes after each adjustment of the whip. Since several adjustments are 
normally required, the process can be quite time-consuming. 

Adjustments are particularly time-consuming if one is receiving 
one of the Omega stations. Each of these stations transmits only about 
one second out of 10, so the long integration time problem is compounded. 
At the time of the experimental work for this paper, there were four 
Omega stations in operation (there are planned ultimately to be eight). 
Each station operates on a time-shared basis on two frequencies, 10.2 
and 13.6 kHz. They transmit an unmodulated waveform in a prescribed 
sequence for a prescribed length of time (e.g., 0.9 sec, 1.1 sec, etc.) 
to facilitate identification. Still in the planning stage for the sys- 


tem are the addition of 11.33 kHz by all stations on the time-shared 
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basis, as well as transmissions of 12.00 and 12.75 kHz (not on a time- 
shared basis) by some of the stations. In order to obtain phase and 
amplitude information on an individual station, it is necessary to 

clamp the AGC and blank the receiver during the time the desired station 
is not transmitting. In the laboratory at the Naval Postgraduate School, 
this was accomplished by a Tracor Omega Gating Unit. This unit may be 
used with four tracking receivers simultaneously, and a modification for 
the receiver is provided which allows switching between normal receiver 
operation and Omega operation. The unit proved highly satisfactory and 
much valuable data was obtained using Omega stations. 

It was useful in the experimental work to have an auxiliary VLF 
receiver available which could conveniently be tuned from station to 
station for quick audio identification. The receiver used was a model 
AN/SRR-11, and a similar type would no doubt be useful for auxiliary 


purposes in the shipboard VLF relative navigation system. 
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5. Frequency Standards 

A good frequency standard must be used at both the phase-stabilized 
VLF station and the ship using the VLF relative navigation system; hence 
it is a vital part of the system. Its use in the System aboard ship is 
shown in Figure 4-1. ' 

There are two basic types of precise frequency standards: the quartz 
crystal oscillator and the atomic frequency standard. The former is 
characterized by a predictable frequency drift as the crystal ages, 
while the latter has either an extremely small drift or none at all. 

The past twenty years have seen great achievements in the field 
of precise frequency control; today, time difference or frequency may 
be measured more precisely than any other physical quantity. Using an 
atomic hydrogen maser, a frequency stability (RMS deviation from the 
mean) of 2 parts in 101^, over a one-day averaging period, may now be 
obtained [30]. One problem which has arisen in the growth of technol- 
ogy of frequency standards is a lack of standardization of the defini- 
tions of certain terms. Two of these are short-term stability and 
long-term stability. The former refers to average frequency deviation 
measured over a time interval sufficiently short that long-term effects 
may be neglected. Since the averaging period is not standardized, it 
must be specified; it varies from a fraction of a second to perhaps one 
minute. Long-term stability refers to slow changes in average frequency, 
and is usually specified for crystal oscillators by stating fractional 
frequency deviation per day, e.g., 5 parts in 1011 per day. Since 
atomic frequency standards have little or no systematic drift, specifi- 
cations for them will often be given in terms of standard deviation or 


RMS deviation from the mean over a specified time interval such as one 
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Quartz Crystal Oscillators 

There are many factors which affect the frequency of quartz crystal 
oscillators, which must be taken into account or compensated for. In 
particular, frequency Stability is a function of: (a) temperature; (b) 
time (aging); (c) stress, vibration, and acceleration; (d) drive level, 
and (e) nuclear effects. These effects are discussed by Gerber and 
Sykes [29]. 

Ambient temperature changes account for the principal frequency 
change effects in most crystal units. Current practice in precision 
oscillators is the use of a double proportional oven to control the 
temperature of the crystal unit. This provides a method of removing 
heat from around the crystal at the same rate it is generated. A sim- 
pler, but less desirable, method uses a thermostatically-controlled 
oven. In this method, temperature varies between a lower and an upper 
limit, whereas it is maintained constant when using proportional con- 
Ero b- 

The effect of time (aging) in a quartz crystal oscillator is easily 
observed by the variation in frequency deviation from day to day. Accord- 
ing to Gerber and Sykes (29), aging of thickness shear crystals is caused 
by four factors: (a) temperature gradient effects lasting several min- 
utes to several hours after a thermal disturbance, (b) stress relief 
effects, (c) change of mass effects, and (d) structural changes due to 
imperfections in the crystal lattice. In all crystal oscillators there 
seems to exist an initial stabilization period, normally lasting a few 
weeks, in which the oscillator has not yet attained its ultimate fre- 


quency stability. 
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The effects on frequency stability of stress, vibration, accelera- 
tion, etc., may be reduced by careful design of mounting structures. 
Austin [24] has made an experimental study of the effects of vibration, 
shock, and swinging motion on frequency stability of several commercial 
crystal oscillator units. He concluded that stability of these units 
was degraded negligibly under conditions similar to those which might be 
encountered in a shipboard environment. 

Another factor affecting crystal stability is drive level, or ampli- 
tude of vibration. A typical 2.5-MHz crystal unit might exhibit a change 
in frequency deviation of from 1 part in 10? to about 8 parts in 107 with 
a change in drive level from 3 to 1000 microwatts [29]. For this reason, 
automatic gain control is used in precision crystal oscillators. 

Nuclear radiation also has an effect on the frequency of crystal 
oscillators; the extent of this effect is discussed in the literature 
[13], [15], апа [26]. 

Figure 5-1, from Gerber and Sykes [29], shows a schematic diagram 
of a modern precision crystal oscillator. The double oven is shown, and 
proportional control would be used. To maintain a constant drive level, 
a large amount of negative feedback would be used in the amplifier, in 
addition to the AGC. 

The Sulzer Model D2.5 is a crystal oscillator of the type shown in 
Figure 5-1 and was in use at the Naval Postgraduate School for the ex- 
perimental work of this paper. Its drift rate was determined to be 
about 2.5 parts in 1011 per day. Other units of this general type often 
have drift rates of about 1 part in 1011 per day, while similar units 
used to stabilize transmissions of some U.S. Navy VLF transmitters have 


average daily drift rates as low as a few parts in 1013 [29]. This 
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Figure 5-1. Precision Crystal Oscillator 
general type of oscillator should be well suited for shipboard use 
in a VLF relative navigation system. 
Atomic Frequency Standards 
It is interesting at this point to note the difference between 
a primary: and a secondary frequency standard. The former is a self- 
calibrating device which is intrinsically accurate. The latter, 
however, requires calibration dc, 0e" a known precise frequency 
before it can be considered a standard. A quartz crystal oscilla- 
tor may thus be classified only as a secondary frequency Standard. 
There are several types of atomic frequency standards, including the 
hydrogen maser, the cesium beam standard, and the rubidium vapor 
standard. The latter, considered a secondary frequency standard 
since its frequency initially must be compared with some other 
Standard, was used at the Naval Postgraduate School as an adjunct 
to the preparation of this paper. The hydrogen maser and the cesium 


beam standard are both primary frequency standards. All three of 


these atomic standards are now commercially available. 
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Atomic standards in general derive their precise frequencies from 
frequencies corresponding to the energy involved in the transition of 
electrons between energy states in particular types of atoms, such as 
hydrogen, cesium 133, and rubidium 87. The cesium and rubidium standards 
operate on the same basic principles; both have found wide use in VLF 
applications as well as other uses. Space does not permit detailed des- 
criptions, but since the rubidium standard was used in some of the ex- 
perimental work for this paper, its basic operation will be described. 

The heart of the rubidium standard is a high-quality quartz oscilla- 
tor, whose frequency is stabilized by a passive resonance cell filled 
with rubidium and an inert buffer gas. The operation of the resonance 
cell is based on an optical pumping principle. The rubidium vapor in 
the cell is illuminated by filtered light from a rubidium lamp. A 
microwave signal is applied to a microwave cavity surrounding the reson- 
ance cell. This signal is obtained by multiplying the frequency of the 
quartz oscillator. When the frequency corresponds to atomic resonance, 
energy transitions are induced in the rubidium atoms and a photo detec- 
tor observes these transitions through changes in the amount of light 

/ 

absorbed as a function of applied microwave frequency. A servo loop 
connects the detector output and oscillator, locking the oscillator to 
a frequency corresponding to the center of the resonance line. Since 
the resonance frequency is influenced to some extent by the pressure of 
the buffer gas in the resonance cell, this pressure must be properly 
adjusted when the standard is manufactured. The cell in then sealed 
and the unit is ready to serve as a frequency standard. 

Rubidium vapor frequency standards have been designed for U.S. Army 


tactical applications and for guided missiles and aircraft. They should 


also be readily adaptable to a shipboard environment. State-of-the-art 
characteristics for these standards, according to McCoubrey Lan). in- 
clude an RMS deviation from the mean of 5 parts in 1012 for a one-day 
averaging period and a systematic drift of less than 3 X inne per month. 
Cesium beam frequency standards are presently being used as the 
United States Frequency Standard (maintained by the National Bureau of 
Standards, Boulder, Colorado) and to stabilize the VLF transmitters of 
the Omega navigation system stations. While the short-term stability 
of cesium beam standards is not as good as that of rubidium standards, 
their principal advantages are excellent long-term stability and high 
intrinsic reproducibility. McCoubrey [30] gives the state-of-the-art 


-13 per day (RMS deviation from the mean). 


long-term stability as 2 X 10 
Stations of the Omega navigation system effect an intercomparison of 
their cesium beam standards periodically (currently every four hours), 
and are thus able to control their stability to better than 1 X 10711 
per day [32]. 

The Varian Model V-4700 is the rubidfum vapor standard that was 
used at the Naval Postgraduate School for the authors' experimental 
work. Over a period of about four months, this unit exhibited an RMS 


deviation from the mean of 1.36 X 10711, using an averaging period of 


one day. 
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6. Antennas 

As mentioned in Section 3, the long path signal at times can intro- 
duce serious error in the VLF relative navigation system. For this 
reason, it is desirable to use some means of discriminating against this 
component of a received VLF signal. The cardioid antenna system, which 
has been widely used in radio direction finding and radio compass appli- 
cations, provides a ready solution to the problem. The basic theory of 
this system has been described by Terman [3] and others, and will be re- 
iterated and elaborated upon in this section. The loop antenna, basic 
to a cardioid set-up, will first be discussed, then the cardioid system 
itself, and finally the results of measurements on antennas at the Naval 
Postgraduate School will be presented. 

The Loop Antenna 

The loop antenna, in various forms, has been used for many years. 
One of its earliest applications was reception of Stations in the broad- 
cast band when this Service was in its infancy. Loops used for VLF 
work, as well as most other applications, have a diameter which is small 
compared to a wavelength. A loop may be circular, rectangular, tri- 
angular, diamond-shaped, or some other similar configuration. Kraus [2] 
States that far fields of loop antennas of different shapes are the same 
if the areas enclosed by the loops are the same. 

A vertical loop antenna produces a polar radiation pattern in the 
familiar figure-eight shape, as shown below in Figure 6-1. The xy-plane 
is parallel to the surface of the earth. Maximum response, of course, 
occurs when the plane of the loop is oriented parallel to the direction 


of arrival of the electromagnetic wave. 
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Figure 6-1. Loop Radiation Pattern 

Appendix III describes the theory of operation of the loop antenna. 
Equations are derived there for the relative response and phase shift as 
the antenna is rotated from 0° clockwise to 360°. Figures 6-2 and 6-3 
show plots of these equations for 0° to 180°. The frequency is taken 
to be 20 kHz, and the loop is assumed square, with the length of a side 
equal to 1.635 m. (This length was chosen because it was the measured 
length of a loop available for experimental work.) The loop is assumed 
parallel to the y-axis of Figure 6-1 and direction of arrival is from 
right to left along the x-axis as shown. The angle of rotation & is 
taken as the angle between the normal to the loop and the positive x-axis 
(this is shown in Figure III-2), so that Ф = 0° and 180° corresponds to 
zero response, while Q = 90° and 270° are the points of maximum response. 

The graph of Figure 6-2 was arbitrarily terminated at -35 db for 
convenience, since output theoretically goes to - о0о db at & = 0° and 
1809. For rotation from 1809 to 3609, the amplitude pattern is exactly 
the same as for the first 1809 of rotation (this is obvious from in- 
spection of Figure 6-1). 

For the phase pattern, the second 1809 was not plotted in order to 


show greater detail. At àQ = 180° and Q = 360°, there is a discon- 
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Figure. 6-2, Loop Relative Output 
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Figure.6-3. Loop Relative Phase Shift 


49 


tinuity of 180°. The pattern for © = 180° to & = 360° looks the same 
as Figure 6-3, except it appears inverted, i.e., at ® = 180°, Q — 
-180°; at & = -270°, 9 = -(180° + 3.4 x Toy deg); and at Q = 360°, 
(Q—- -1809. (The symbol —» is used to imply the existence of the dis- 
continuities at & = 180° and & = 360°.) The 180° phase reversal on 
the back side of the loop antenna is significant, for one must be care- 
ful not to rotate the loop antenna to the wrong side of the null if 
phase meaSurements are important. 
The Cardioid Antenna System 

A cardioid antenna system consists of a loop antenna and a whip 
(or "sense") antenna whose outputs are appropriately combined to pro- 
duce a cardioid-shaped radiation pattern. Figure 6-4 is a simplified 
schematic of such a system. Both the loop diameter and the height of 
the whip are small with respect to a wavelength. Some provisions for 
tuning, which are not shown, would normally be installed at the base of 
the loop and at the base of the whip. There would also be provisions 


for varying the amplitude of the loop signal or the whip signal, or both, 
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Figure 6-4. Cardioid System Schemat ic 
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as well as provisions for varying the whip phase a small amount above and 
below the 90° phase shift shown. These provisions are necessary to allow 
adjustment to obtain the cardioid pattern. The adjustment consists of 
making the whip antenna component of the VLF signal equal in magnitude 
and phase to the loop contribution, for the case in which the plane of 
the loop is oriented toward the distant transmitting station. 

Appendix IV contains the equations necessary to plot the phase and 
magnitude of the cardioid antenna system as the loop is rotated about a 
vertical axis. The necessary phasor diagrams are shown, and Figure IV-3 
shows the phasor diagram for proper cardioid adjustment, with the loop 
oriented toward the station. Proper operation of the cardioid system is 
independent of the location of the whip with respect to the loop, so long 
as it is not an appreciable fraction of a wavelength away. 

The radiation pattern of a properly-adjusted cardioid antenna sys- 
tem is sketched in Figure 6-5. Maximum response occurs with the plane 


of the loop antenna oriented parallel to the x-axis. 


4 


DIRECTION OF 
SIGNAL ARRIVAL 


Figure 6-5. Cardioid Radiation Pattern 
Figures 6-6 and 6-7 are plots of the radiation pattern in Car- 
tesian coordinates, using equations derived in Appendix IV. The plot 
is for a frequency of 20 kHz and a square loop 1.635 m on a side.  Ro- 


tation angle & is the angle between the normal to the loop and the 
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positive x-axis. The maximum response is 6 db, which merely indicates 
that the voltage delivered to the receiver is twice that from the loop 
or whip alone. The minimum response, with the loop pointed away from 
the direction of signal arrival, is approximately -63 db. In the forward 
direction, 6 db beamwidth is 1809, providing the desirable characteristic 
of a broad front lobe. 

It should be noted at this point that -63 db is not the minimum 
response obtainable on the "back" side of a cardioid antenna system. 
The ideal minimum response approaches -oO0db when the system is adjusted 
for this. The adjustment required for this entails adjusting the whip 
phase for minimum system response, which corresponds to setting the angle 
of E. to *(99; rather than -@ (see Figures IV-4 and IV-7). The reason 
for the present approach to the adjustment problem is that for a weak 
station it is nearly impossible accurately to adjust the whip for a good 
null while the antenna is pointed in the reverse direction; this has been 
verified in the laboratory. 

Phase shift of the cardioid antenna system is plotted in Figures 
6-8 and 6-9, again using equations derived in Appendix IV. It should 
be noted that the Y-axis scales of the two figures are different. For 
0°<&< 180°, the phase shifts no more than 0.3 x 107% deg as the loop 
rotates. At a frequency of 20 kHz, 360° phase shift corresponds to 50 
psec, so 0.3 x 10-2 deg would correspond to less than 107? usec; since 
normal phase variations of propagation produce changes of 1 or 2 Msec 
on a phase plot of a received VLF station, these small changes would not 
be noticeable. 

An interesting aspect of Figure 6-9 is the phase change of -90° in 


the vicinity of Q = -90°. At 20 kHz, this would correspond to a phase 
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change of 12.5 psec. It should also be noted that, while there may appear 
to be a discontinuity at @ = 0° and & = “71809, due to the difference 
in Y-axis scales on Figures 6-8 and 6-9 this is not in fact the case. 

Since it is proposed that a cardioid antenna system be used aboard 
ship in the VLF relative navigation system, there has been some concern 
that the sharp phase change in the vicinity of the cardioid null might 
cause difficulties. It is believed by the authors that this will cause 
no problems for the following reasons: (a) the relative broadness of the 
null, and (b) experimental evidence that the 90° phase shift in a typical 
cardioid system does not occur so abruptly as the theoretical curve of 
Figure 6-9 indicates. The latter will be discussed later in this sec- 
tion. The former is seen by examining Figure 6-7: if the loop is rota- 
ted 15° either side of the null, cardioid response is still down nearly 
35 db from that of the forward direction. For these reasons, it is felt 
that no special measures would need to be taken aboard ship to keep the 
antenna aligned with the direction of the transmitting station, other 
than a simple servo system driven by the ship's gyro. 

Experimental Work 

The Tracor Model 611 Cardioid Unit was used at the Naval Postgraduate 
School. This unit has two inputs: one directly from a loop antenna, and 
the other from a whip antenna. Both antennas must be tuned as desired in 
the VLF band and matched to 50-ohm coaxial cable. The cardioid unit pro- 
vides adjustable phase shift and amplification so that the cardioid pat- 
tern may be obtained at any frequency between about 9 kHz and 30 kHz. 
The output of the unit is switchable between cardioid, loop only, and 
whip only. Although the cardioid unit is specified to be usable with 


only particular types of Tracor antennas, it was found to be usable with 
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an 18-foot whip antenna that happened to be on hand and a tuning network 
designed for 20 kHz. The complete cardioid system manufactured by Tra- 
cor comes with a remotely-tuned rotatable loop antenna and a whip 
antenna with fixed tuning. This system was used to obtain data pre- 
sented below. 

Radiation patterns for the cardioid system were obtained using WVL 
at 20 kHz. The length of a side of the square loop was 1.635 m, so data 
is directly comparable with the theoretical patterns of Figures 6-6, 6-7, 
6-8, and 6-9. The system was adjusted by setting whip amplitude and phase 
equal to that of the loop with the loop oriented toward Ft. Collins, 
Colorado (the location of WWVL). It would have been preferable to set 
up a nearby low-powered (milli-watt or micro-watt) phase-stabilized 
transmitter to assure steady signal levels while the antenna was rotated, 
but this was not practicable. 

The amplitude data obtained is shown in Figure 6-10. Azimuth of 09 
corresponds approximately to orientation toward the transmitting sta- 
tion. The numbers in parentheses below the abscissa correspond to angle 
of rotation Q. Due to difficulty in reading the rotator unit, relative 
azimuth accuracy is probably +30. Relative amplitude is accurate to 
better than 21 db. The two major sources of error here are variations 
in the received carrier level of WVL and calibration of the receiver 
carrier level recorder output. 

It is interesting to compare Figure 6-10 with Figures 6-6 and 6-7. 
Considering that real antennas never perform quite so well as ideal ones, 
the similarity is striking. Overall, the null is somewhat broader than 
the theoretical; for example, the output is 35 db down at about 22° either 


side of the null. The 6 db beamwidth is 1819, which is the same as the 
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theoretical value, within the experimental error. The two points appear- 
ing within the null at amplitudes of -32.5 db and -35.4 db indicate that 
the cardioid system was misadjusted, with whip amplitude slightly great- 
er than loop amplitude. A careful examination of Figure IV-7 shows that 
this would also cause phase excursion to be slightly greater than ~A/4. 
This was verified experimentally, as seen in Figure 6-11. 

Figure 6-11 shows phase shift data obtained as the loop was rotated. 
The ^f = 0 Msec line is an arbitrary reference level. Phase data is 
probably accurate to Íl or 2 sec, considering accuracy of the receiver 
and propagation changes that could have occurred while data was being 
taken. The smaller figure is probably acceptable, since data was taken 
during the daytime when the ionosphere was stable, over a period of 
about three hours. Angle of rotation & is again indicated in parenthe- 
ses below the abscissa. The most notable difference between this curve 
and the theoretical data of Figures 6-8 and 6-9 is the broader width of 
the 90° (A/4) phase change in the vicinity of the null, as noted ear- 
lier. The rubidium vapor frequency standard was used as the local phase 
reference for this data. 

There are various problems associated with using a cardioid antenna 
system in a shipboard environment. One is the problem of re-radiators 
located near the loop-whip system. This presents a problem in a receiv- 
ing system because of distortion of the antenna pattern, principally by 
resonant lengths of metal near the antenna. Unfortunately, it was not 
possible to investigate this problem at the Naval Postgraduate School. 
It is, however, the opinion of the authors that this problem would not 
seriously degrade the cardioid antenna's performance in the VLF relative 


navigation system. 
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Cardioid Phase 


Figure 6-11, 


Another problem is that of properly adjusting the cardioid unit. 
Whip phase and amplitude must be carefully matched to that of the loop 
at each frequency used, in order to obtain a true cardioid pattern. This 
may be difficult for distant stations whose long and short path signals 
are comparable. Careful attention should be given these adjustments, 
since depth of the null depends directly on how well the system is ad- 


justed. 
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7. Determination of Signal Strengths 

In this section values will be calculated for the relative strengths 
of signals received via the long and short great circle paths as a func- 
tion of the time of day and day of the year. This, along with the phase 
of both signals, is necessary in order to calculate the phase of the re- 
sultant signal. 

The relative strengths of the two signals can be found at any point 
on earth if the total attenuation on each path is known. The relative 
strength of the two signals is then just the difference of the short and 
long path attenuations. Upon reexamination of the electric field strength 
equation in Section 2 it will be seen that an assumption is being made 
here that all terms except the last two in this equation are identical 
for both propagation paths. It seems only reasonable to assume that most 
VLF stations have a symmetrical radiation pattern and that the amount of 
excitation in two opposite directions will be equal. 

Attenuation rates at VLF, mentioned in Section 2, vary from about 
one to five db/Mm depending on the time of day and propagation path con- 
ditions. References 5, 8, 9, 18, and 28 contain many of the curves of 
attenuation rates versus frequency which have been derived from experi- 
mental data. These curves will be used extensively here. 

The procedure used to obtain expected attenuation rates at any par- 
ticular frequency consists simply of sifting through the data present in 
these references and extracting the attenuation rates given for different 
conditions at this frequency. From this information estimates of the 
attenuation rates for all possible:/ path conditions and times of day can 
be made. 


As an example, at 13.6 kHz the following values are given in the 
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references stated. 


CX (db/Mm) 


1 O 


Ax = 1.0 to 1.5 


A07 1.0 


A&K= 1.0 


conditions 
land, day, east to west 
sea, day, east to west 


day 


sea, night, east to west 


day, sea, west to east 


night 


sea, day, east to west 
sea, day, north to south 
sea, day, west to east 
night, sea 


sea, night, west to east, 
15-20 kHz 


sea to land change, 
15-20 kHz 


night to day change, 
15-20 kHz 


W-E to E-W change, 
15-20 kHz 


reference 
Taylor [9] Ре 933 
Тау1ог [9] ре. 1933 


Chapman and Macario [5] 
pg. 930 


Taylor and Lange [8] 
pg. 609 


Taylor [9] pg. 1933 


Chapman and Macario [5] 
pg. 9230 


Watt and Croghan [18] 
Watt and Croghan [18] 
Watt and Croghan (18) 
Watt and Croghan [18] 


Davies (28) pg. 425 
Davies [28] pg. 425 
Davies [28] pg. 425 


Davies [28] pg. 425 


` 


From these attenuation rates (and changes in attenuation rates) and 


other discussions of attenuation rates such as those by Crombie [78] and 


Westfall [21] reasonable estimates of attenuation rates for 13.6 kHz can 


be made. These will be in the form of one basic attenuation rate plus 


additional factors which depend on the direction of propagation, amount 


of daylight on the path and amount of land on the path. 


The values chosen are: 
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e, - 1.0 db/Mm basic rate (night, sea, W-E) 
Ахра = 1.0 db/Mm night to day change 
& X. 4 - 0.8 db/Mm sea to land change (average land) 
AG. s = 0.4 db/Mn W-E to N-S direction change 
А Сару = 1.55 db/Mm W-E to E-W direction change 


The change in attenuation rate as a function of direction can be 


approximated by: 


AX (g) = 1.26 @ 2 


1809 
where Ø is the direction of propagation with respect to 090° true 
(0 < D 4 1800). For north to south propagation @ = 90° and for east 
to west propagation @ = 180°. 
The attenuation rate over any path can be calculated from the fol- 


lowing equation: 


ОС = HK, + PAK, 4 + PLA, 1 AXC) 


where Pp = fraction of the path which is in daylight 


fraction of the path which is over land. 


M 

A typical plot of attenuation for a 24 hour period is shown in 
Figure 7-1. The path here is Trinidad to Monterey (short path); @ = 
исо, Pr = 0.6 апа Pp is obtained from the computer program described in 
Appendix II. The total attenuation is found by multiplying & by the 
path length obtained from the geodesic program (Appendix I). The cal- 
culated attenuation rates vary from 2.48 db/Mm for the entire path in 
darkness to 3.48 db/Mm for the all-daylight path. 

Figure 7-2 shows a comparison of the long and short path attenua- 
tion from Trinidad to Monterey. The long path here is considered to be 


in the opposite direction at the receiver and transmitter from the short 
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UT 
Calculated attenuation of 13.6 kHz signal on 
Trinidad to Monterey short path for 1 Oct. 1966 
Path length ы 6658.85 km | 
Nine 16.220 
Max.=23.37 ар 
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Figure 7-2. Calculated attenuation of 13.6 kHz signal on 
Trinidad to Monvepey Lone ng. shorts pachsei or 
4 Oct. 1966 | 
Short path length z 6658.85 km 
Long path length = 33391.48 km 
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path direction of propagation. 

It should be noted that the attenuation rates used here for land 
paths are for average land. Large areas of very high or very low conduc- 
tivity land would require alteration of these values. Similar care must 
be exercised for propagation over polar ice caps, where the conductivity 
is lower than that of sea water. 

One factor not mentioned previously is the loss due to discontin- 
uities of the earth conductivity along the path (La n in the electric 
field strength equation). Westfall [21] mentions this factor briefly 
as having a magnitude of three to six db at times. Because of the lack 
of information on this factor it has not been taken into account here. 

The method shown in this section could be used at any other fre- 
quency of interest in the VLF range to obtain similar attenuation rates. 
Some complicating factors may be encountered, however, which could make 
the calculations more difficult. For example, at 16 kHz it was found that 
the night to day difference over land is greater than the night to day 
difference over sea. 

Comparison of Experimental and Calculated Values 

Figures 7-3 and 7-4 show two more curves of calculated attenuation 
for the Forestport and Haiku Omega stations and Figures 7-5 through 7-7 
show the observed relative carrier levels of these three stations as ob- 
tained from a Tracor Model 599G VLF phase tracking receiver. No measure- 
ments were obtained on any long path signals or on the Aldra short path 
signal because the carrier levels were below the level of control of the 
receiver AGC, giving a reading of zero on the relative carrier level out- 
put. 


Since no attempt was made to calculate the signal strength of the 
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Figure 7-3. Calculated attenuation of 13.6 kHz signal on ~. 
Forestport to Monterey short path for i Oct. 1966 
Path length = 4011.74 km 
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UT 
Figure 7-4, Calculated attenuation of 13.6 kHz signal on 
Haiku to Monterey short path for 1 Oct. 1966 
Path length = 3856.50 km 
Min. = 4.16 db 
Max. = 8,02 db 
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Figure 7-5. Observed relative carrier level of Trinidad 
13.6 kHz ot Monterey, 30 Aug. to 2 Sept. 1966. 
10 db/ in. vertical scale 
night average - day average: = 4.04 db 
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Observed relative carrier level of Forestport 
beet kHz signal st Monterey, 5 Sept. to 8 Sept. 
1966. 

10 db/in. vertical scale 

night average = day average = 2. 27 ар 
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various stations and this information is not readily obtainable from the 
Tracor receivers, the only quantities that can be compared meaningfully 
are the difference between the carrier levels for the entire path in 
darkness and in daylight and the day-night attenuation difference. While 
the attenuation decreases at night and the carrier level increases during 
the same period the difference between the day and night values should 

be equal since they are complementary quantities. (This is easily seen 
in the electric field strength equation of Section 2 where c«d/109 is 


the attenuation and Es is the carrier level.) 
3 


n 


A comparison of the calculated and observed values of the night-day 


difference is shown below. 


Station Calculated difference Observed difference 
(a) (b) 
Trinidad 6.65 db 4.04 6.4 db 
Forest port 4.01 db Qd. 7. b 
Haiku 3.86 db 0.3 4.5 db 


The values listed under (a) were observed with the receiver servo 
time constant set at 150 seconds and those under (b) with a 15 second 
time constant. The lower differences observed with the 150 second time 
constant, corresponding to Figures 7-5 and 7-7 are apparently because 
of the erratic nature of nighttime phase. 

Calibration of the signal level scale used on the recorder was very 
difficult, leading to some possible inaccuracies in these measurements. 
(Calibration was complicated by the short duty cycle of the transmitting 
stations; 10%.) Thus while these measurements do provide some measure of 
confidence in the calculated attenuation rates, it is believed that more 
comprehensive measurements of the carrier levels of these signals would 


be necessary if an accurate check of these attenuation rates was desired. 
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8. Calculation of VLF Phase Velocities 

Figures 2-8 through 2-11 show several curves of calculated phase 
velocities for VLF signals for different conditions which could be used 
directly to obtain the desired phase velocities. However, it would be 
more convenient and accurate to make use of an equation for phase velocity 
which takes into account the frequency and conditions of the ionosphere 


and the earth. Such an equation is available in the Omega System book 


[31], page 3-20: 








V h; у 2 
P u= 1 0936 ee | opm - (0 :) ——— 
А : Vr ФО тту 
Неге Vp = phase velocity of the VLF signal 
Vg ^ free space velocity of the wave 


= 3x 108 m/sec 


0.36 = a constant chosen to give a good fit in the 10-20 kHz 
frequency range 


h., = ionospheric height for which the conductivity (g) = 1076 


mho s/m 
a = radius of the earth = 6378.155 km 
Фе = phase shift upon reflection at the ground 
Ui - phase shift upon reflection at the ionosphere 
n = mode number (= 1 for primary mode) 
f = frequency 
Some values are given for the phase shifts that can be expected; 
(p. - -1609 for daytime propagation, -1659 at night, -165° for daytime 
propagation to the east, and “155 гог daytime propagation to the west. 
Assuming that the -1609 phase shift is for north-south propagation, which 
seems reasonable, the values for nighttime propagation to the east and 


west would then be -170° and -160° respectively. The ground phase shift 
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is zero for sea water and the values given for land are -0.09 radians 
during the day and -0.13 radians at night, for a ground conductivity of 
1072 mhos/m. 

Using the given ionospheric heights of 68 km during the day and 86 
km at night (derived from an exponential conductivity profile and taking 
the height for which 7 = 107? mhos/m) the phase velocity curves of Fig- 
ure 8-1 were calculated. A comparison of these curves with similar 
curves in Section 2 (Figures 2-9, 2-10, 2-11) shows a very close corres- 
pondence. Several experimental phase velocities from Figure 2-14 are 
also plotted here for further comparison. 


The phase velocity at any point along a path is: 


2 2ТЕ 
Ури). = чиле 
where © is the phase in radians at the point x on the path. The total 
phase along the path is then: 
: d 
S mij VE 
O 


Vp (x 


If the phase velocity is constant this becomes: 


с оаа 
Ө у 


P 
Using a method similar to that employed in determining the path 


attenuations of Section 7, an average phase velocity can be obtained 


for each path. It can be calculated by the following equation: 


Ыз Aa (B) + PPAV р + РУ 1, 
where Урь = basic phase velocity; phase velocity for night, land, 
W-E from Figure 8-1. 
А Урас) = increase in phase velocity due to direction of propagation. 
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Pp = proportion of path in daylight. 


АУ Р-р = increase in phase velocity for daylight propagation. 


Р 


S proportion of path over sea. 


âVp1 increase in phase velocity for propagation over sea. 


-S 
The increase in phase velocity with direction can be obtained from 


the following equation: 


AVpa P) = AVpa 1509 


where @ is the direction of propagation, measured as in Section 7. From 
Figure 8-1, the differences between the phase velocities for E-W and W-E 
propagation are different for daytime and nighttime propagation. &У ьа 
for daytime can be approximated by the average daytime differences in 
phase velocity with direction and the nighttime value can be obtained 


similarly. This can be expressed by the following equations: 
— ë 
А үа = %(VpsDEW ` Ypspwe * Ypipew ~ Ypipwe? daytime 


= 1 _ ° ° 


Here V = phase velocity for sea, day, E-W 


psDEW 


V = phase velocity for land, night, W-E and so on. 


рІММЕ 


Using reasoning similar to the above it is found that two values 


for 2, must be used (one for E-W propagation and one for W-E) and 


N-D 


one value for АУ 1-5 can be used. These are given by the following 


equations: 


V 


ae x: E E 
4 Von-p y "CV sDWE VosNWE ў Vp1DWE pINWE? "AP 


sah E = 
A V n-p = @%pspew ` VpsNEW * Ypipew ~ Vpinzw > E-W 


V V 


m ali m = 
AVp1-s B 4) s DWE Vp 1DWE i psDEW p 1 DEW 
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VDSNWE ^ 'plNWE * VpsNEW ^ Yp1NEw? 

The range of the phase velocity changes being averaged to obtain 
^ Vp1-s is about 4 x 107? V Yo" The range of phase velocity changes 
being averaged for АУ N-D and ^ Vpà is smaller, indicating that the 
error in & V51.s is the largest that would be encountered. 

The curve shown in Figure 8-2 is a typical 24 hour plot of phase 
delay calculated by the above method. This is for the Trinidad to Mon- 
{егеу short path with the same parameters as used in Section 7 (Ø = 


1459, Bess 0.4, Py from computer program of Appendix II). The 13.6 kHz 


phase velocities used to calculate this curve are listed below. 


VospWwE = 1.0000213V, 
VospEw = 1000427 3V, 
Vosnwe = 9-9974367V, 
Vosnew = 0° 9976842V, 
Vpipwg = 0-99981997, 
Voipew = 1-0002154v, 
Voinwe = 0-9972606V, 
Voinew = 0+9974986V, 


Comparison of Experimental and Calculated Phase Delays 
Figures 8-3 through 8-6 show both calculated and observed 24 hour 
phase delay curves. The average of the observed daily curves is shown 
shifted by 1.5 inches (15 or 30/Asec) for convenience and ease of com- 
parison with the calculated curves. The total phase delay of the signal 
from transmitter to receiver is not easily obtained experimentally; thus 
only the day-to-day changes in the phase delay are compared here. 


Figures 8-5 and 8-6 show a very close correspondence between calcu- 
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lated and observed curves. It is interesting to note that there appears 
to be a time lag in the observed phase delay during periods of rapid 
change, which could be due to the long time constant used in the receiver 
(150 sec). The lag also seems more pronounced on stations with lower 
Signal strengths at the receiver. 

Figure 8-4 shows differences in the calculated and observed curves 
which seem to point towards some error in the calculations. The ampli- 
tude and shape of the observed diurnal shift are believed to be approx- 
imately correct. 

Figure 8-3 is a good example of a diurnal shift that is not a trap- 
ezium. Propagation over a polar path,as is the case here, yields 
diurnal shifts similar to those shown in Figure 8-7, and are markedly 
different during different times of the year. The aisi shift on polar 
paths is discussed more thoroughly by Blackband (12). In Figure 8-3 the 
poor correlation of the calculated and observed curves may be partially 
due to the weak signal strength of the Aldra station at Monterey. The 
phase-tracking receiver did not seem to "lock on" to the received phase 
until a large change in phase delay occurred, such as encountered during 
sunset along the path. This can be seen during the initial portion of 
the readings, started at 1900 on 23 September. 

Figure 8-8 shows a comparison of an experimental phase velocity 
curve, obtained by Thompson et al [14] in Sydney, Australia at 16 kHz, 
with a calculated curve. The data used for the calculated curve is 
approximate (position of the receiver, proportion of land on the path, 
etc.) and in view of this the close correspondence of these two curves 


is remarkable. 
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Figure 8-5. Forestport to Monterey phase delay, short path 
13.6 kHz, 6-9 October, 1966 


81 


Phase — 









Eg ee 
i E sh 
iet zo | Seen 
es ‚4 1 
hase -$ +- рф 9 ё + 
, m ae iiu. š 
a t 


: L. 
Pre tee 


al H 














! i 
| 2 ЦЕ ТЕП 
т Г 
ar 










ay & Hom 

pred лес ae 
a ts 9 uS 
d TTE Я a E IE de T , [ 
























ae 
A 
d e p l M 
Tm SEQUENCE AES 


FEE rH -HÁ БЕЕН ыг It Eu =m A XA ШЫ 


ШЕКЕН И 
б к 





PE EP Ee F 
о. т 
у ч aa 


Бо + 
H MARA 








сз 














t ГЕТ te SI sS i [sj ыш ү еа n 


Г тыр Ime pee 5L CL ds 
» em [сл Bt 





















a SEHE +N Ñ HE 


= p 
== THR ЕЕЕ = Hee ЕН: НҮН 


eS 
See Tee Te us 


. == ut o E 
sss t UM Е AES 
Е 


He zm Be "dept jn 
| F TH 
= = cr us Es Gene ee iri p Hr — = EE сЕ са = ас Site A? 
Se UE ES eee ee eee reer EH ius ЕНЕВ 


ЕЕЕ ЕЕС СЕ ЕЕ 
Eo e EM шш | | |o c oom 


Ex I ERO UE Ei PIE E 
SEE E s shi sns EE 


n a2 à 6 8 10 312. i14 16 18 20 22D 
UT 


EIE s 








mcm m THH 


Figure 8- 6. Trinidad to Monterey phase delay, short path 
13.6 kHz, 11-14 October, 1966 
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9. Calculation of Resultant Signal Phase 

As discussed in Section 3, the amplitude and phase of a signal 
consisting of long and short path components can be obtained by the 
addition of two phasors. In the navigation system proposed in this 
paper the phase of the resultant signal at a receiver is of primary 
interest. The amplitude of the signal is of interest only insofar as 
it 1s sufficient to allow phase tracking. 

The phase of the resultant signal can be found from the phase of 
the long and short path signals and their relative amplitudes or atten- 
uations. First take the magnitude of the phasor Ë In Figure 3-3 to be 
unity, or zero db. Then the magnitude of f, is the relative attenuation 
of the short path signal with respect to the long path signal. It can 
reasonably be assumed that when lE I< 0.1 the phase of E, is not appre- 
ciably affected by E, (for €, = 0.147909, Ф = 5.79 in Figure 3-3). 

The relative attenuation for this case is given by the following 
equation: 


E 
A = 20 log tot = - 20 db. 


Thus the long path signal must be taken into account whenever the 
strengths of the two signals are within approximately 20 db of each 
other. 

To obtain the phase of the resultant signal, the phase delay in 
microseconds, calculated in Section 8, must be converted to an angle. 
This is easily done by converting to wavelengths, where one wavelength 
equals 1/f seconds (at 20 kHz, 1/f » 50/4sec). One wavelength corres- 
ponds to 3609 and 8 is then the difference between the phase delays on 
the long and short paths. The angle Ø is then given by the following 


equation: 
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This angle could also be found by adding graphically the phasors Ë) and 
Fn. 

It must be remembered here that the total phase delay of the re- 
sultant signal is equal to Ø plus the phase delay of the short path sig- 
nal. 

Figures 9-1 through 9-5 show the calculated curves of long and short 
path attenuation, and short path and resultant phase at Sydney, Australia, 
and Wellington, New Zealand. The transmitter station used was GBR at 16 
kHz and the day was assumed to be 23 September 1966. 

In Figure 9-2 the largest difference between the short path phase 
and the resultant phase is slightly less than one microsecond and the 
difference in attenuations at that time is 20 db. This verifies the 
earlier assumption that the resultant phase differs from the short path 
phase (by an appreciable amount) only when the strength of the long path 
Signal is within 20 db of the short path signal. 

Figure 9-4 shows a case where the resultant phase is considerably 
different from the short path phase. The largest difference here is 12.1 
asec (а 0300) or about 70 electrical degrees, and does not occur at the 
time when the long path signal is strongest (long path signal is 0.6 db 
weaker at 0300, 1.8 db stronger at 0500). From about 0800 to 1800 there 
1S approximately 20 db of separation between the two signals and the re- 
sultant phase is very close to the short path phase as expected. 

In Figure 9-5 the latitude of Wellington was changed by four min- 
utes, lengthening the short path by about seven km and shortening the 


long path by the same amount. The total phase delay on the short path 
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Flgure 9-2. Short path and resultant phase of GBR at Sydney. 
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Figure 9-4. Short path and resultant phase of GBR at Wellington 
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Figure 9-5. Short path and resultant phase of GBR at 
Wellington. (Short path 7 km longer) 
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is increased slightly as a result of this, but the diurnal shift remains 
the same. The resultant phase, however, is considerably different, show- 
ing that a change in position has a larger effect on the resultant phase 
than on either the long or m ar phase alone. 

No experimental data has been found with which these calculated 
curves could be directly compared; however, Thompson et al. us) and 
Crombie [7] have published articles which discuss this phenomenon and 
include some experimental evidence taken at Sydney and Wellington. 

It should be noted that the resultant phase curves as shown here 
could not be observed by using a loop antenna. The 180° phase reversal 
with a reversal of the direction of the incoming signal would cause 
some different resultant phase to appear. This could be taken into 
account in the calculated curves, however, by adding 180° to the phase 
of the long path signal before adding the phasors. 

The concept of phasors as discussed here provides another means of 
describing the phenomenon of cycle-slipping. Figure 9-6 shows the am- 
plitude and phase of some station over a 24-hour period as represented 


by a phasor trace. 





Figure 9-6. Phasor representation of a 24-hour 
phase and amplitude curve 
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This figure is drawn assuming a trapezoidal diurnal pattern of both 
phase delay and amplitude, both being greater at night. Thus À repre- 
sents the daytime phase and amplitude, B represents the nighttime phase 
and amplitude, ý is the diurnal phase shift, and |B| - |A| is the night- 
time increase in amplitude. 

Cycle-slipping has occurred if the curve traced by the phasor en- 
closes the origin, as shown by the dotted line in Figure 9-6. This 
could be due to long path interference, loss of signal (in which case 
the magnitude of the vector is zero), or several other causes. The 
case for which the diurnal shift is greater than one сус1е (21) should 
not be confused with cycle-slipping as shown here. Figure 9-7 shows one 
possible curve of this type and it can be seen here how easily cycle- 


slipping could occur. 





Figure 9-7. Phasor representation of a 24 hour phase and 
amplitude curve with a large diurnal phase shift 


A rapid change in phase, as sometimes experienced during sunset or 
sunrise, could cause the carrier level in the receiver to be greatly re- 
duced for a time, particularly if a long servo time constant is used. 
This could result in a phasor trace similar to the one shown in Figure 


9-8, with an increased probability of cycle-slipping. This rapid phase 


9:3 


change and resultant cycle-slipping 1s typical of north-south propagation 


paths. B 


Figure 9-8. Phasor representation of loss of carrier 
level due to rapid phase change 
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10. Conclusions 

Attenuation rates and relative signal strengths of long and short 
path signals can be determined easily. The accuracy of the method used 
here has not been determined, however, because of a lack of dependable 
experimental data. 

It has also been demonstrated that the phase of a signal can be 
predicted for a typical 24 hour period, both with and without long path 
interference. Experimental data has shown the short path phase calcula- 
tions to be very accurate for some cases when three-to five-day averages 
of observed phase are considered. Lack of experimental data prevented a 
check of the calculated phase in the presence of long path interference; 
however, it was shown that when this type of signal is present the phase 
is highly sensitive to a change in receiver position. 

Additional experimental data, possibly from some location other 
than Monterey, California, would be highly desirable and would probably 
lead to some modifications of the parameters used in the calculations 
presented here. It should be emphasized that these calculations are a 
first iteration in an attempt:to predict, at any point on earth, the 
phase of a VLF signal. 

A cardioid antenna system should be suitable for use aboard ship 
in a VLF relative navigation system. A simple servo system driven by 
the ship's gyro, and provided at intervals with the true bearing to the 
desired VLF station, should be adequate to keep the cardioid pattern 
oriented with sufficient accuracy. With proper adjustment of the car- 
dioid system, interference from the long path signal may be eliminated. 
In ease of a location at which the long path signal is always equal to 


or greater in strength than that of the short path, however, it may be 


25 


impossible to adjust the cardioid system properly, and hence it could 


not provide the desired amount of long path signal attenuation. 


¥ 
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APPENDIX I 
GEODESICS AND LONG PATHS 
A.  Geodesics 
The great circle distance between two points on a sphere is easily 
calculated; however, determination of the length of a geodesic (analo- 
gous to a great circle) on an ellipsoid is more difficult. 
Kirkland [20] states very concisely the simplest method of obtaining 
the length of a geodesic. 
Given two points (6,, А у), (Po; r,) find their respective para- 
metric latitudes, 8 у апа B 2 as shown in Figure 1-1. 
Ø = geodetic latitude 
= parametric latitude 


а. = equatorial radius 
b = polar radius 











tan @ = - -3z 
dy b 
Z=a cos 8 
y =b sing 
dz a 
= =- —— tan 
dy b Ê 
_ b _ 
tang = — tan @ RNC tan ó 
sec?43 - 1 + c? tan? @ Figure I-1. Geodetic and 


Parametric latitudes 


1 
cos б = 
sec 


cos 
cos? ø + c2 sin? ø 


sing = Vi - соѕ24 = c San 


cos? (b + с? sin? ø 
The arc distance between the two points (x) is then found from: 


соѕ х = зіп 3 1 зіп 8 + соз @ ү соз 8 5 соз (Ду - А 2) 
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and the great circle distance is D=ax. 
Now using the correction factors derived by Lambert [1] the length 
of the geodesic is: 


d = D - AP - BQ 


(a - b) (x - sin x) 


where P = i 8 ESSO 
que Dr Ges) 


4(1 - cos x) 


> 
I 


o ^ (in(8, * sin (8 5? 


Bo (sin f$ i - sin @ ))2 

The accuracy of this method is, according to Lambert, better than 
+ 0.05 km but breaks down near the antipodes (i.e. for x >175°). This 
method was checked against the examples in an article by Rainsford [4] 
with the largest error being about 0.85 km. (It should be noted that 
Rainsford specified the position of the points to within 0.0001 second.) 
B. Long Path Distances 

The method used in this paper to obtain the long path distance be- 
tween ы points is to find the length of the perimeter of the ellipse 
on which the two points lie (analogous to the great circle on a sphere) 
and subtract the geodesic length from this distance. 

The major axis of this ellipse is the equatorial radius of the 


ellipsoid and the minor axis can be found as follows: 


8 


parametric latitude 


x = arc distance between the two points 
cos @% = cos A - cos x cos C 


Sin x sin C 
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зіп З > - cos x sin /$ | 
sin X COS З 





EQUATOR 
Figure I-2. Ellipse of geodesic 


| sin($5 - cos x sin($j 
cos § = sin (х cos 6 = cos 1 1 - [ a: | 
1 B sin x cos /3 ү 


Here gamma is the parametric latitude of the northernmost point on the 


ellipse. From Figure I-3: 


2 = а cos% 
22 у? 
+ ac xu 
а^ b b 
2 
т = b2 m > z2 = b2(1 - соз? y ) 
a 
f- 


comici Ж 
а. 
Figure I-3. Parametric and 
geocentric latitudes 











The minor axis of the ellipse of interest is then: 


r = 4p? + y? = -/ъ? + (а? - b^) cos? y 


The length of the perimeter of the ellipse described by a and r is 


found by evaluating the elliptic integral: 


/2 
ЕСК) = A - к^ sin2x dx 
о 


where 


k = 4/1 - (— 2 
a 
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The perimeter is then: 
P=4 aE (k) 

and the long path distance is: 
d, =P - d 

where d is the geodesic length. 


The FORTRAN 60 program used to compute these distances is shown on 


the following pages along with sample results. 
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APPENDIX II 
PROPORTION OF DAYLIGHT AND DARKNESS ON PATHS OF PROPAGATION 

The difference in the height of the D-layer of the ionosphere from 
day to night is the largest single cause of variations in the phase vel- 
ocities of VLF waves and causes variations in the attenuation rates of 
a magnitude equal to those caused by other factors. Because of this the 
percentage of daylight and darkness on each path of interest must be 
calculated as a function of time. 

This problem has been solved by Brady and Crombie [22]. However, 
since the following solution was developed independently, it will be out- 
lined here. 

The following is based on the premise that the percentage of day- 
light on a path can be found if four points are known: (1) transmitter 
location, (2) receiver location, (3) and (4) intersection of the sunrise 
and sunset lines with the great circle path between (1) and (2). (A 
spherical earth of unspecified radius is assumed.) 

The great circle path can be described by the longitude of one of 
the points where it crosses the equator CÀ) and its declination at 
the equator (O«,). Тһе sunrise-sunset line forms another great circle 
that can be described in the same way. The intersection of these two 
great circles on the sunrise side of the earth will determine the coor- 
dinates of point (3) mentioned previously. The geometry of the situation 
is shown in Figure IL-l. 


Here: 


^2 


15t - 90 degrees 


(C = K, + 90 degrees 
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where 







sunrise line 


t = GMT ircle path 


declination of the sun 


x 





equator 
2. 
Figure II-1. Point of sunrise on 
great circle path 


Using spherical trigonometry the coordinates of P СА}; ^ 3 can 


easily be found: 


cos4s(&, - OX j) 
a = tan”! | tank A, z Z 


cos%( x) + Ko) 


+ 


sink(€, -w» 1 
vans Е jt À шл. 


sin’s(&K 1 * 6€») 
tan c = tan a cos(180° - > ,) - - tan a cos CX» 
sin É 3 = sin a sin %> 


The point at which the sunset line crosses the great circle path 


is then: 


^, 
B, - -f; 


At a height of 70 km the point P will be further along the great 


+ О 
À * 180 


circle path in a direction away from the sun. Figure II-2 shows the 
are distance by which the sunrise and sunset lines are shifted. 
Here: 


h = ionospheric height 
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2 
li 


screening height 


ye] 
li 


earth radius 


m 


sun's rays 
Figure II-2.  Sunrise/sunset at ionospheric height 
The screening height, mentioned by Brady and Crombie [22], is the 
height below which the components of the sun's rays that most affect 
lonization cannot penetrate. This, in effect, changes the radius of 
the earth by an amount S. From Figure II-2 then: 
Razas 
+ 


cdi 2 


The point at which this new sunrise line crosses the great circle 
path (P4) is shown in Figure II-3. The coordinates of this point can 
be found by spherical trigonometry: 


cos*s( X 4 - x 2) 
cos%(% 1 + O 2) 


tanl(a -* d) 


tams( À> ч À 1) 


i sinz(X1- X33 
сап А, = Мр р атару 


ёап5(а =- а) 


; gi 
sins(a + d) 
з} wl 2 з1п%(а - d) 


: » sin Y 
sindAd = P y 
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j t 
Figure II-3. Revised position of sunrise point 


From Figure II-4: 


cos P = sin x cos( À 3 - A 1) 


cot Ad cos 63 - sin; cos A 
sino 


cot Ө 


^3 


Аз + Ө 


sin @ 3 =” IT 4 3 cosd + соз з зїпда cos О 





Figure II-4. Sunrise position coordinates 
The coordinates of the point of sunset (i> Ai) can be found in 


a similar manner from Figures II-5 and II-6: 


Au = 3 £ 1809 


s - 83 


^d = -Ad 
ee cot Ad cos By + sin@y cos (o 


sin Ø 
А» = А» + е 
sinó, = sin, cos Ad - cos ц sin Ad cosf 
yos A EE 
e 





^u 





Figure II-5. Revised position of sunset point 





40^ -8 
i S PoLE 


Figure II-6. Sunset position coordinates 
It should be noted here that if it is assumed that É and Аа 
are negative quantities, then @' will also be negative and the resulting 


longitude correction will be in the opposite direction from the sunrise 
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longitude correction. 

Once the coordinates of the four required points are determined 
the amount of sunlight can easily be found with the aid of Figure II-7. 
This shows the four possible situations that can exist on the great 
circle path of propagation. Points (1) and (2) are the transmitter and 
receiver positions, (3) is the sunrise point, and (4) is the sunset 
point. The amount of daylight on the short path for case A is the dis- 
tance from (3) to (2); the amount of daylight on the long path for case 


A is the distance from (4) to (2), and so on. 


ILSR (4) 55 


e (CN p 
S. D (I) 


(4) 55 (3)5R 


case A case B case С case D 





Figure II-7. Possible positions of sunrise and sunset 
on great circle path 


The actual program is shown on pages 116 through 120. The geodesic 
program was converted to a subroutine for convenience here. In the pro- 
gram the declination of the sun as a function of the time of the year 
is approximated by two sine waves of different periods; one being used 
from the lst to the 80th and the 267th to the 365th days, and the other 
being used from the 8lst to the 266th day. Data from the Nautical Al- 
manac [177] was used as a guide in setting up this approximation. 

The program also includes a correction factor for the Greenwich 
Hour Angle, which was also found necessary from data in the Nautical 
Almanac. This correction factor has a maximum value of 4.1 degrees 


(this changes ^» in Figure II-1) and is again approximated by sine 
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waves of various periods and amplitudes (see statements 311, 511, 711, 
and 812 in the program). 

Geocentric latitudes, obtained from the geodesic program, are used 
so that no error is involved due to the ellipticity of the earth. The 
radius of the earth is not even used in this program, since the distances 
between points are expressed as arc distances and the amount of daylight 
is expressed as a percentage of the long or short path. 

The only inputs required are the time of day and the day of the 
year. The latitude and longitude of the two points of interest and 
the description of the great circle path through these points are ob- 
tained from the geodesic subroutine. 

In the paper by Brady and Crombie [22] three basic simplifications 
were made which could result in some error. These were the assumptions 
that: 

1) the sun is a point source, 

2) atmospheric refraction can be neglected, 

3) the earth is spherical. 

The third assumption does not introduce any error in the program devel- 
oped here because geocentric latitudes were used. The other two could, 
according to Brady and Crombie, result in errors of about 16 minutes 
and 34 minutes of arc respectively in the quantity Ad (Figure II-3). 
These errors, although small, could cause considerable error on high 
latitude great circle paths. 

This program is not valid for the case in which the highest latitude 
on the path is greater than the quantity (909 - E This causes a por- 
tion of the path to be continually in daylight or darkness, (e.g., polar 


paths during the winter or summer) and the method used here breaks down. 
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APPENDIX III 
LOOP ANTENNA CALCULATIONS 
Consider Figure III-1 below. This is a schematic representation 
of a loop antenna showing a voltage generator in each vertical side of 
the loop to indicate the voltage induced by the passing wave. The z- 
axis is taken as the vertical axis and the xy plane the horizontal plane 


with respect to the earth. The voltage Ep apPearing across the tuning 





Figure III-1. Loop Antenna 


capacitor C is the signal produced by the antenna and is found by phasor 
addition to be re - E, = Еу. The symbol ~ here denotes a phasor quan- 
tity, while a letter without this indicates magnitude only. 

In order to find the effect on the antenna output of various di- 
rections of arrival of the signal, assume the antenna to be rotatable 
about the z-axis. The angle & is defined as the angle between the 
normal to the loop and the direction of arrival of the signal, as shown 


in Figure III-2. 


t 
LOOP 


DIRECTION OF 
M 
x ARRIVAL 


Figure III-2. Loop Geometry 


TZ 


The angle varies from 0° to 360° and taken to increase positively as the 
loop rotates clockwise. The bar shown on the loop in Figure III-2 is 
for reference and corresponds to side 2 of the loop as seen in Figure 
III-1. 

The phasor diagram of the loop voltages is given in Figure III-3. 
Ihe angle ( represents the phase of the antenna output signal, where 
(9 = 0 when the loop is oriented broadside to the direction of arrival 


of the wavefront. 





Figure III-3. Phasor Diagram 
The angle 8 in Figure III-3 represents the difference in phase be- 
tween the voltage induced in side 2 and that in side 1 of the loop. It 
1S determined by the distance d between the two vertical sides, and is 
given by 


8 = —— sin& radians, where 


A 


free-space wavelength of the received signal 
(same units as d). 


Since E, = E), the phasor diagram forms an isosceles triangle and there- 


fore Ф = 590. The phase angle Y is then 


p =- T 


sinQ radians. 
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The resultant Ер is also found readily from the phasor diagram. 
Using the Law of Cosines, one has 
Её = ET + ES - 2E1E2 cos Ө 


= 282 = 281 cos Ө 

L 2 

= 2E] (1 - cos 0) 
Ер = Еу V2 (1 - соз Ө) 


= Еу ү2 (02 sin — 


мез? sin Q. 
It is convenient to normalize by setting E, = 1. One can then express 


the magnitude of the relative output Eout in db below ERmax: 


T ER 
out 7 ?9 logio T 


2 sin Q E 
a "е ‘Tsing,’ po Yo Фо 


sin 
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APPENDIX IV 
CARDIOID ANTENNA CALCULATIONS 
As in Appendix III, assume the direction of travel of the incoming 
signal fixed and the loop antenna rotatable, with the angle © defined 
as before; see Figure IV-1. In addition, the sense antenna is placed 
at some arbitrary position, forming the angle (X' as shown. The distance 
d. of the sense antenna from the loop midpoint is non-critical, so long 


S 


as it 1s short with respect to a wavelength. 






SENSE 


DIRECTION OF 
SIGNAL TRAVEL 





| 
Figure IV-1.  Cardioid Geometry. 

Because of the 1809 phase reversal of the loop antenna when 
a = 0° and 180°, it is convenient to analyze the cardioid system in 
two parts, with 09 4 (L4 1809 and with 18094 A < 360° (or equivalently, 
-180°< 064 09). As in Appendix III, the loop phasor resultant Bp is 
given by E» - Ei and phase is 6/2. This is illustrated in the phasor 
diagram of Figure IV-2. In this diagram, both Ө апа W are taken to be 
negative angles. 

For proper adjustment of the cardioid pattern, the magnitude and 
phase of the signal from the sense antenna must be made equal, respec- 
tively, to the magnitude and phase of the signal from the loop antenna 


when the loop is positioned for maximum response (Q = 90°). This 
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Figure IV-2. Loop Phasor Diagram 
condition is shown in Figure IV-3, where E, is the magnitude (after the 
necesSary amplification and attenuation) of the sense antenna Signal. 
The angle (0. is the largest (negative) value of tf. Once set, the 
magnitude and phase of Е, of course, remain constant as the loop is 


rotated. 


$ 





Figure IV-3. Cardioid Adjustment 
The phasor diagram for some value of O. between 0° and 180°, other 
than O&= 90°, is shown below in Figure IV-4. The overall signal out- 
put of the cardioid unit is given by E,/y? (note that ^Y is negative and 
its magnitude less than that of Ф)- This is obtained by the phasor 


relationship Bs = Ер + E,, accomplished in the cardioid unit which mixes 


the loop and sense antenna signals. 
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Figure IV-4.  Cardioid Phasor Diagram 


As in Appendix III, the angles Y, s 8, and 8, may now be ex- 


pressed in terms of the angle of loop orientation Q , as follows: 





Ө = 29 зіп О; 3600 (degrees) 
A 
-d sing 
= —— —Ao @v x 1800 
A 
ONG = m 
о шах 
-d 





0 = Qu - ` x 1809 


The objective is to find “f and E in terms of Q. At this point, 
it is convenient to derive an expression for Ep as a function of Es 
Ihis is done as follows, referring to Figure IV-2 and using the Law of 
Cosines: 


2 2 2 
E Ey + ЕЛ - 2E, E, 


R 


cos 8 


= (1 - cos Ө), since By = Е). 


oa Ep = E, V2 sin c 


Ф 
II 


sin 8, since @ is small. 


ў Ер. 
^ TED e v2. 
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ERmax TM Es 
ЕС ОЛКЕ EQ ` 
Ер _ @ 

Es Ө, 


Figure IV-5 is an expanded view of the phasor diagram of Figure 
IV-4, drawn on an exaggerated scale to make angular relationships 


clearer. In this diagram and the following equations, the fact that 





Figure IV-5. Angular Relationships 


and Ф„ аге negative is taken into account. In the triangle of Fig- 
ure IV-5, the two sides Ep and Eg and the angle between them are known, 
so that J may be found by geometry. One can then determine Ү= ф- 5, 
which is the phase angle of E. Solving the triangle, 


o 1 ES - E 
e л с” E + tan? | 








2 Е. + Ер 
O = 
cot ( КЁ, ) 
2 

S 

F 
- 909 - 909 «Las + tan? — 

1 + 
8. 


cot (909 - — 
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where x = 


Since and Ç can be found for any given value of Q, Ү is uniquely 

determined. It is a negative angle and (p, «'Y < 0° for 00° <Q < 1800. 
Determining the relative amplitude response for these values of (X 

is very straightforward. Refer to Figure IV-5 and use the Law of Cosines 


as follows: 














ЕД = ES + E? - 2EgE, cos | 1809 + CY, - Ф) | 
= E + Е2 + 2ЕрЕ5 cos (49, - 9) 
ЕЕЕ eos (Po -YD 
EL I TR mz S )2 + 2 E cos( Q5 - Y ) 





= 41 + x2 + 2xcos (Qao -ф) 


One may now normalize by setting E, = 1 and taking 20logi oF, to obtain 


relative output in db, with zero db occurring at Q = 0° and 180°: 
E, = 20 logipV1 + x2 + 2x cos((9, - Ф) db. 


The procedures for obtaining cardioid phase and amplitude for loop 
positions 1809 4 (€ 4 3609 are quite similar to the methods given above. 
Figures IV-6 and IV-7 below are self-explanatory and are analogous to 
Figures IV-2 and IV-4, respectively, discussed above. 

Angles are Similar to those determined before, the principal dif- 
ference being that 180° is now added to (Q due to the difference in loop 
position. Note that all of the following angles are now positive, while 


Y remains negative: 
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Figure IV-6. Loop Phasors (1809 < Х < 3609) 





Figure IV-7. Cardioid Phasors (180° < Q < 360°) 


8 = 


Ф = 


Ө 


о 


As before, S 


Ў = 


- 3 sin x 360° 
A 
1809 - c sin(& x 180° 


d 
—— x 1809 
A 


E. 2590. 
A 


is found by trigonometry as follows: 


T о Е Е 
909 - CON RETE tan | | CE 
2 Es + Ep 
+ - 1809 
cot est LES 


Чо + 1 | х-1 Чо + Ф 
о ГЕШ ee ЕНЕ ЦРНЕ с I — Kue n 
180- - 5 + tan =I tan 5 Р 


Y| = 180° - [ç ф- 1809) + x] 
ЖП = | 1809 - [c9 - 1809) + 5] 
= -180° +  - 1809 + X 

Y= P+ 3 - 360°, 


Relative amplitude is found as before using the Law of Cosines, and 


referring to Figure IV-/7: 


[т] 
№ 
II 


2 
2 = E2 « Eg - 2EgEg cos [1809 - (&, +¥)| 


B = М1 + х2 + 2х соз (фо + Фф). 





Normalizing and expressing the result in db, one obtains: 


E, - 20 logjg V 1 * x? * 2x cos(Q, + Ҹ) а. 
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